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Abstract: In arid areas, ecological degradation aroused by over-exploitation of fresh water, expansion of 
artificial oasis and shrinkage of natural oasis, has drawn attention of many scholars and officials. The water 
and ecological footprints can be used to quantitatively evaluate the water consumption of social-economic 
activities and their influence on the eco-environments. In addition, increase of the water footprint indicates 
the expansion of artificial oasis, and the influence on the natural oasis could be reflected by the variation 
of the ecological footprint. This study was conducted to answer a scientific question that what is the 
quantitative relationship between the expansion of the artificial oasis and the degradation of the natural 
oasis in the arid environments of Xinjiang, China. Thus, based on the social-economic data, water 
consumption data and meteorological data during 2001—2015, we calculated the water and ecological 
footprints to express the human-related pressure exerted on the water resources and arid environments in 
Xinjiang (including 14 prefectures and cities), and explore the relationship between the water and ecological 
footprints and its mechanism by using the coupling analysis and Granger causality test. The results show 
that both the water and ecological footprints of Xinjiang increased significantly during 2001—2015, and the 
increasing rate of the ecological footprint was much faster than that of the water footprint. The coupling 
degree between the water and ecological footprints was relatively high at the temporal scale and varied at 
the spatial scale. Among the 14 prefectures and cities examined in Xinjiang, the greater social-economic 
development (such as in Karamay and Urumqi) was associated with the lower coupling degree between the 
two footprints. Increases in the water footprint will cause the ecological footprint to increase, such that a 1- 
unit increase in the consumption of water resources would lead to 2—3 units of ecological degradation. The 
quantitative relationship between the increases of the water and ecological footprints, together with the 
intensities of water consumption both in the natural and artificial oases of Tarim River Basin, have approved 
the fact that the formation and expansion of 1 unit of the artificial oasis would bring about the degradation 
of 2 units of the natural oasis. These conclusions not only provide a technical basis for sustainable 
development in Xinjiang, but also offer a theoretical guide and scientific information that could be used in 
similar arid areas around the world. 
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1 Introduction 


Water is the source of life, the key element of production activity and the foundation of eco- 
environment (Chen, 2012). Rational utilization and optimal allocation of water resources are the 
primary prerequisites for achieving sustainable development in arid areas (Gao et al., 2018). The 
scarcity, irreplaceability and imbalance of water resources in spatial and temporal distributions 
fundamentally restrict the speed of social-economic development in the oases of arid areas, which 
are inseparable from freshwater resources (Jiao et al., 2013; Gao et al., 2018). The amount of water 
consumption represents the developmental scale of regional oases (Liu et al., 2010; Chen et al., 
2019). Prosperity or ruin of oases can be attributed to the abundance or lack of water, and many 
ancient civilizations on the Silk Road (such as Loulan and Niya) failed due to lack of water 
resources and over-exploitation of natural resources (Qin et al., 2012; Zhang, 2016). In addition, 
there are plenty of contemporary examples worldwide that regional environmental degradation in 
arid areas was caused by excessive water use, such as Aral Sea Basin, Heihe River Basin and Tarim 
River Basin (e.g., Liu et al., 2010; Zhang, 2016; Gao et al., 2018; Chen et al., 2019). Therefore, 
unreasonable exploitation and utilization of water resources are regarded as the main causes of 
environmental degradation in arid areas. 

Ecological footprint analysis, proposed by Rees (1992) and improved by Wackernagel and Yount 
(1998), is an effective method to evaluate whether the current relationship between human beings 
and natural resources is sustainable or not (Wackernagel et al., 1999; Xu et al., 2003). The 
ecological footprint analysis approach reflects the most important advance in the quantitative 
evaluation of sustainable development. This approach attempts to define sustainable development 
using measurable indicators and to monitor the depletion of natural resources caused by humans. 
The ecological footprint, which contains the components of built-up land, carbon, cropland, fishing 
grounds, forest products and grazing land, can quantify all the resources consumed by a specific 
population for production activities and the total ecological area needed to eliminate all the waste 
generated by the population. These objectives of the ecological footprint analysis are achieved by 
establishing a comprehensive accounting system that compares human ecological need and 
ecological carrying capacity of eco-environment, and explores the way of implementing sustainable 
development (Rees, 1992; Wackernagel and Yount, 1998). If the ecological carrying capacity is 
larger than the ecological footprint, the ecosystem is sustainable, otherwise the ecosystem is 
unsustainable and humans have to change the current mode of development to make it sustainable. 
The ecological footprint has been widely applied on various scales, from global to national, 
community to family, and commercial enterprise to individual activity, for its advantages of vivid 
concept, abundant connotation and comprehensive theory. 

Similar to the ecological footprint, the water footprint was developed to assess the water 
consumption by human activities (Ridha and Ben, 2018). The water footprint is defined as the 
amount of water needed for any product and service consumed by human beings at the country, 
regional or individual scale over a certain period of time (Hoekstra et al., 2011; Veettil and Mishra, 
2016). The products and services here include foods and daily necessities for human life, domestic 
water and environmental water (Bosire et al., 2015). The water footprint analysis involves 
systematic evaluation of the impacts of human consumption activities on water resources in terms 
of the quantity, type and efficiency of utilization. It also clarifies the demand for water resources 
and enables rational evaluation of the utilization of water resources in both economic production 
and daily life. Therefore, the water footprint is mainly applied to evaluate the security, carrying 
capacity and sustainable utilization of water resources (Zhang and Li, 2012; Veettil and Mishra, 
2016). 

Recently, a series of integrated footprint approaches were put forward to assess the 
environmental impacts of production and consumption activities. For examples, Galli et al. (2012) 
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developed a concept of "footprint family" in the One Planet Economy Network: Europe 
(OPEN:EU) project and selected the water footprint, the ecological footprint and the carbon 
footprint to define the concept. Ewing et al. (2012) proposed a multi-regional input-output 
framework for the integrated calculation of the water and ecological footprints. Fang et al. (2014) 
conducted a theoretical exploration by combining the water footprint, the ecological footprint and 
the carbon footprint. To sum up, integrated footprints (Hertwich and Peters, 2009; Rushforth et al., 
2013; Sun et al., 2013; Hao et al., 2015) were widely used in assessing the environmental 
sustainability (Fang and Duan, 2015) and ecosystem services (Zhang and Zhang, 2013), evaluating 
the safety of drinking water (Qi and Chang, 2012) and exploring variations in the water, ecological 
and carbon footprints at different scales. However, most previous studies focused on the conceptual, 
theoretical and physical applications of the footprints, researches about the connections and 
mechanisms within the integrated footprints were scarce. 

Indeed, the high consumption of water resources cannot only promote the social-economic 
development, but also damage the natural environments, as reflected by the expansion of the 
ecological footprint, especially in ecologically fragile region such as Xinjiang Uygur Autonomous 
Region, China. The total water footprint of Xinjiang ranked the 7" in China and the per capita water 
footprint in Xinjiang was as high as 1400 m*/a, more than twice of the national average level in 
China (Ge et al., 2011). Synchronously, the total ecological footprint of Xinjiang ranked the 25" in 
China, and the per capita ecological footprint was 2.91 ghm? (global hectare), which was slightly 
higher than the national average level in China (2.86 ghm?) and the global average level (2.82 ghm?) 
according to the statistics from the Ecological Footprint Network 
(https://data.footprintnetwork.org/#/). There is a consensus that over-exploitation of water resources 
by human activities will inevitably result in the degradation of environments (Yang et al., 2006; 
Salameh, 2008; Zhang et al., 2017), which could lead to the high water and ecological footprints. 

In arid areas, oasis is the main area for the development of human production and consumption 
activities. The artificial oasis is the core area that human lives in, and the natural oasis spreads 
around outside the artificial oasis to prevent it suffering from desertification. The areas of both 
natural oasis and artificial oasis are approximate in Xinjiang, and in terms of modern land use 
dynamics, the natural oasis has been continuously transferred to the artificial oasis (Wu et al., 
2014). Taking the Tarim River Basin as an example, with the expansion of the artificial oasis, the 
water amount consumed by human activities has risen sharply, thus resulting in a gradual decrease 
in the amount of water entering the terminal lakes. Further exacerbating the situation, the lower 
reaches of the Tarim River has dried up for almost 30 a from 1972 to 2000, severely affecting the 
local environment and making some areas no longer habitable; as a result, a housing transfer 
phenomenon happens. Thus, the dynamics of the artificial and natural oases have drawn attentions 
of many researchers. Suitable scale of oasis (Hu et al., 2007; Ling et al., 2012, 2013) and 
appropriate area proportion of the natural and artificial oases (Li et al., 2008; Guo et al., 2017) in 
arid areas have been studied. Yet, the mechanism of transformation from the natural oasis to the 
artificial oasis and the influence of expansion of the artificial oasis on the natural oasis have not 
been revealed clearly. The increase of the water footprint indicates the expansion of the artificial 
oasis (Chen et al., 2019), and its influence on the natural oasis can be reflected by the variation of 
the ecological footprint. Understanding the development of the water and ecological footprints is 
essential to explore the mechanism between the dynamics of the natural and artificial oases. 

The objectives of this research are to analyze the relationship between the water and ecological 
footprints, assess the quantitative influence of water over-exploitation on ecological degradation 
and reveal the mechanism and quantitative relations between the expansion of the artificial oasis 
and the degradation of the natural oasis. To achieve these goals, we calculated the water and 
ecological footprints of Xinjiang from 2001 to 2015, analyzed the connection between them at the 
spatial and temporal scales and explored the Granger causality relationship between the two 
indicators. This study not only provides a technical basis for sustainable development in Xinjiang, 
but also offers a theoretical guide and scientific information that can be used in similar arid areas 
around the world. 


2 Study area and data sources 


2.1 Study area 


Xinjiang Uygur Autonomous Region is the hinterland of the Eurasian continent and has a total area 
of 1.6610° km?, accounting for almost one-sixth of China's land area. A total of 14 prefectures and 
cities are included in Xinjiang: Aksu Prefecture (abbreviated as Aksu), Altay Prefecture (Altay), 
Bayingol Mongolian Autonomous Prefecture (Bayingol), Bortala Mongolian Autonomous 
Prefecture (Bortala), Changji Hui Autonomous Prefecture (Changji), Hami Prefecture (Hami), 
Hotan Prefecture (Hotan), Kashgar Prefecture (Kashgar), Karamay City (Karamay), Kizilsu Kirgiz 
Autonomous Prefecture (Kizilsu), Tacheng Prefecture (Tacheng), Turpan City (Turpan), Urumqi 
City (Urumqi), and Ili Kazak Autonomous Prefecture (Ili). Xinjiang locates in a temperate arid 
region with the mean annual precipitation of only 146 mm and the mean annual potential 
evaporation of 1600-2300 mm. The annual mean temperature in Xinjiang is 7.9°C, and it is relative 
high in plains and basins and low in mountainous areas (Li et al., 2013). The total oasis area is 
0.1610° km’, occupying 9.6% of the total area in Xinjiang (Wu et al., 2014). The average annual 
river runoff in Xinjiang is 8.79x10!° m? and the available water resources is 8.35x10!° m? (Deng, 
2010). 

Xinjiang is the most prominent region where the scarcity of water resources is severe in China. 
Nationally, the water consumption in 2017 was 5.52x10!° m°, which accounted for 66.1% of the 
total regional available water resources. The regional gross domestic product (GDP) in 2017 was 
1.09x10!? CNY and the added value of primary industry was 0.16 10!? CNY, accounting for 14.7% 
of the regional GDP. Due to the superior in light and heat resources, the agricultural planting 
occupied a relatively large proportion in the national economy of China. The total planting area 
was 6.06*10° hm? in 2017, which accounted for almost 78.0% and 40.0% of the artificial oasis and 
total oasis areas, respectively. Large-scale cultivation of agriculture has sharply risen the 
agricultural water consumption, and the water consumption in primary industry was 5.14*10!° m’, 
accounting for 93.1% of the total water consumption of Xinjiang. 

In 2017, the per capita GDP in Xinjiang was 4.49x104 CNY, which was 75.0% of the national 
level in China. The per capita GDP varied significantly among the 14 prefectures and cities inner 
Xinjiang. For instance, Karamay had the highest per capita GDP, with the value of 15.80x10* CNY; 
Urumqi, the capital of Xinjiang, had the medium-high per capita GDP of 7.78x10* CNY; Kashgar 
and Hotan had the lowest per capita GDP, with the value of approximately 1.07x10* CNY. 
Generally, Urumqi and Karamay had developed industries and low amounts of water consumption, 
whereas Kashgar and Hotan relied mainly on planting industries and the amounts of water 
consumption were relatively high. Therefore, the low-level economic structure has caused a serious 
imbalance in the water consumption structure, further influencing the ecological environments. 


2.2 Data sources 


Social-economic data describing the population, GDP, consumption of major foods by rural 
households, consumption expenditure of urban households, prices for farm products and energy 
consumption of Xinjiang during the period 2001-2015 were obtained from the China Statistical 
Yearbook (National Bureau of Statistics of China, 2002—2016), China Energy Statistical Yearbook 
(National Bureau of Statistics of China, 2002-2016), Xinjiang Statistical Yearbook (Bureau of 
Statistics of Xinjiang Uygur Autonomous Region, 2002-2016) and Xinjiang Production and 
Construction Corps Statistical Yearbook (Bureau of Statistics of Xinjiang Production and 
Construction Corps, 2002—2016). We calculated the water footprint of main crops in Xinjiang: 
wheat, corn, rice, potatoes, beans, cotton, oil-bearing corps, vegetables, melons, apples, pears, 
grapes, jujube, alfalfa, hemp and medicinal materials. The agricultural added values (the increased 
values of agriculture, forestry, animal husbandry and fishery, which are the balance of the current 
value after deducting the intermediate input) at constant price were obtained from the 
aforementioned Xinjiang Statistical Yearbook and Xinjiang Production and Construction Corps 
Statistical Yearbook. The data describing the amount of water resources, regional consumption of 
water resources and agricultural water consumption were obtained from the Xinjiang Water 
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Resources Bulletin (Water Resources Bureau of Xinjiang Uygur Autonomous Region, 2001-2015). 
We divided the oasis into natural oasis and artificial oasis according to the land use types (Fan 
et al., 2002; Wang and Guo, 2018). Specifically, the natural oasis consists of forestland, grassland 
and natural water area; and the artificial oasis is composed of arable land, orchard, built-up land, 
reservoirs and water channels. The land use data with the spatial resolution of 30 m in 2010 were 
obtained from the Geographic Information Monitoring Cloud Platform (http://www.dsac.cn/). The 
climate data (including daily hours of sunshine, average wind speed, relative humidity, 
precipitation, and maximum and minimum temperatures) during the period 2000-2015 for the 
water footprint calculations were obtained from the website of China Meteorological 
Administration (http://data.cma.cn/), which contained 88 meteorological stations in Xinjiang. 


3 Methods 


3.1 Ecological footprint 


Following the approach developed by Wackernagel et al. (2004), we calculated the ecological 
footprint by combining the footprints of renewable resources, fossil-fuel and built-up areas. The 
ecological footprint of renewable resources, also called the biomass footprints (Feng, 2011), was 
used to estimate the area of bio-productive land occupied by humans (e.g., arable land, grassland, 
forestland and water area for fishery). It was determined by using the consumption data of 11 types 
of renewable products, i.e., crops, processed oils, vegetables, fruits, beef, lamb, pork, poultry, dairy 
products, eggs and aquatic products. The footprint of fossil fuels was estimated by using the 
consumption data of coal, fuel oil, crude oil, gasoline, kerosene, diesel oil and natural gas. The 
footprint of built-up areas was quantified as the bio-productive land area occupied by human 
settlements and infrastructures. For a given region, the ecological footprint was calculated using 
the following formula (Wackernagel et al., 1998, 1999, 2004; Yang et al., 2000; Xie et al., 2001): 


EF = Nxef =Nxr;x > (4)=Nxr;x_P/IYN,, (1) 


where EF is the ecological footprint (ghm?); N is the amount of local population; ef is the per capita 
ecological footprint (ghm?/capita); j is the type of land use; r; is the equivalence factor of the j® 
land use type (i.e., the ratio of the global average potential productivity from the j® land use type 
to the global average potential productivity in terms of all bio-productive land use types), which is 
used to standardize different types of bio-productive land; A; is the area of bio-productive land 
calculated by the consumed product i (ghm?); P; is the amount of each product i consumed (t); and 
YN; is the Chinese average consumption of product i, or the carbon uptake capacity in the case of 
the fossil fuel footprints (t). To ensure the continuity of ecological footprint across different years, 
as well as the comparability between our results and the findings in the relevant literature, we 
obtained the constant factors of YN; and 7; from the calculation of Wackernagel and Yount (1998). 


3.2 Water footprint 


The water footprint consists of the crop water footprint (CWF), the animal-related water footprint 
(AWF), the industrial water footprint (IWF) and the domestic water footprint (DWF). The 
calculation of CWF was based on the method described in the Water Footprint Assessment Manual 
(Hoekstra et al., 2011). The CWF was calculated as: 
CWF= SY OWVE sie + WE procgreen) x Y,, (2) 
i=l 
where CWF is the total volume of the crop water footprint (m+); i is the crop category; WF; procsblue 
is the blue water footprint of the crop i (m*/kg); WF; proc,green represents the green water footprint of 
the crop i (m?/kg); and Y; is the yield of the crop i (kg). We calculated the green and blue water 
footprints of crops according to the procedures and descriptions in Hoekstra et al. (2011) and Zhang 
et al. (2018). 
In determining the system boundary of the animal-related water footprint, only the related 
activities substantially contributing to the total water footprint were considered. Therefore, the 


animal-related water footprint primarily consisted of two parts: (1) the total amount (i.e., the 
theoretical amount) of virtual water consumed by the animals throughout the entire life cycle (from 
birth to slaughter); and (2) the total amount of water consumed during the processing of meat 
products after slaughter. 

The total amount of virtual water consumed by the animals throughout the entire life cycle was 
divided into three parts: growth and processing of feed water, drinking water and service water. 
The AWF can be calculated as follows: 


AWF = (VWC,,.g + VWWC ari nk + WW gervice) X P, (3) 


where AWF is the animal-related water footprint (m*); VWCreea, VWCarink and VWCervice represent 
the growth and processing of feed water, drinking water and service water (m*/t), respectively; and 
P is the amount of animal meat products (t). 

The water content in feed materials (i.e., growth and processing of feed water) consisted of two 
parts: the water requirement in mixed feed and the virtual water contained in different feed 
ingredients. The total drinking water supply was equal to the amount of virtual water in the drinking 
water. The amount of service water was equal to the amount of water used in the animal's life cycle 
to clean the housing (e.g., barns or sheep houses), clean the animals and keep the environment 
clean. The amount of water consumed following the slaughter of animals and during the processing 
of meat products was considered to be a part of the industrial production process. 

The processes of formation and consumption of virtual water in the industrial production and 
domestic life were complicated, so we directly assumed the water consumptions by industry and 
domestic life as the IWF and DWF, respectively. 


3.3 Homogeneous indices and coupling degree 


To match the ecological footprint with the water footprint, we conducted dimensionless 
normalization of the two footprints by using the following equation: 
X 


S Re 
i=l 


where X; is the sequence of the water footprint or ecological footprint; Y; is a new sequence (i.e., 
homogeneous indices) of X; after homogenizing; and n is the capacity of sequence Xi. 

According to the concept of capacity coupling in physics and a model of the capacity coupling 
coefficient, we obtained the coupling degree model of multiple systems or factor interactions. For 
analysis, the coupling function between urbanization and the environment can be directly obtained 
and expressed as follows: 


Y. = 


L 


(4) 


C={(u, x uz )/[(u; + Uy)* (u +u), (5) 
where C indicates the coupling degree between the water and ecological footprints; and uı and u2 
represent the dimensionless normalization of the water and ecological footprints, respectively, 
which can be calculated by Equation 4. According to Zhang et al. (2016), the C value is between 0 
(the minimum coupling degree) and 1 (the maximum coupling degree). When the C value equals 
to 0, it implies that the coupling degree is extremely small, the increase of the water footprint is not 
related to that of the ecological footprint, and the system devolves to disorder. The specific division 
of the C values and their implications are shown in Table 1. 


3.4 Granger causality test 


The Granger causality test, firstly proposed by Granger (1969), is an important causal test method 
in economics. Granger (1988) then abstracted the definition of causality, introduced the concept of 
the information set and provided a general definition of causality based on the timing of an event 
occurrence. However, it should be noted that the early definition of causality based on a time series 


Table 1 Coupling degree between the water and ecological footprints 


C value Level of coupling Correlation between the water and ecological footprints 
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0.0-0.5 Low Weak 
0.5-0.8 Middle Moderate 
0.8-1.0 High Strong 


Note: C indicates the coupling degree between the water and ecological footprints. 


is easier to be understood and consequently it is most widely used. The regression of other variables 
(consisting of their own past values), including the hysteresis value of x, can significantly improve 
the prediction of Y. If an estimated or measured entity x significantly improves the prediction of Y, 
x is considered to be the Granger cause of Y. Likewise, Y can be defined as the Granger cause of X, 
which is the sequence of x. The Granger test enabled the establishment of the following model: 


Y, = Soya + Yaa, +c, (6) 
i=l i=l 


where ¢ is the time trend; Y;is the time sequence; y is the element of sequence Y; and n is the lagging 
period. In this study, the water and ecological footprints were in turn fitted to Equation 6 as the left- 
hand side variables. Hurlin and Venet (2001) have assumed that the autoregressive coefficient a; 
and the regression coefficient slope B; are constant, and i belongs to [1, n]. They have also assumed 
that a; is identical for all individuals, whereas the regression coefficient slope B; can have an 
individual dimension. This setting forms the basic framework for Granger causality in a panel data 
context. In Equation 6, c is a constant. 
According to Hurlin and Venet (2001), the homogenous non-causality hypothesis test is required 
to identify the Granger causality between the water and ecological footprints, and it was given by: 
Ho:61 = B2= +: = Ba = 0, (7) 
where Ho (null hypothesis) is the non-existence of causal relationships across n. If the null 
hypothesis is rejected, x is proven to be the Granger causality of Y. In the general case, the F test 
statistic could be computed by using the following Wald test proposed by Hurlin and Venet (2001): 


_ (RSS, —RSS,)/n (8) 

RSS, N —2n-1) ° 

where RSS; is the unrestricted sum of squared residuals computed from Equation 6; RSS2 denotes 
the restricted sum of squared residuals obtained under the null hypothesis; n is the degree of 


freedom (df) of RSS2; and N denotes the total number of observations. Thus, the statistic obeyed 
F-distribution with df (n, N—2n-1). 


3.5 Statistical analyses 


The analyses of correlation and statistical test were processed by software package SPSS 20.0 for 
Windows (SPSS Inc., Armonk, NY, USA), and the Granger causality test was conducted by EViews 
8.0 for Windows (IHS Global Inc., Irvine, CA, USA). 


4 Results 


4.1 Temporal and spatial variations of the water and ecological footprints 


According to the methods stated above, we calculated the water and ecological footprints and 
analyzed their temporal and spatial variations, as shown in Figure 1. In Figure la, the total amount 
of the ecological footprint increased significantly (described by the regression equation: y=0.48x7— 
1902x+1.9x 10°, R?=0.96), from 38.80 10!° ghm? in 2001 to 146.94x10!° ghm? in 2015. The sharp 
increases in the ecological footprints of fossil-fuel land and farmland contributed the most to the 
increase of the total amount of the ecological footprint. Except for the decreasing magnitude of the 
ecological footprint of forestland, the ecological footprints of other components all increased. 

The total amount of the water footprint also increased significantly during the study period, from 
293.12x10® m? in 2001 to 558.47x10° m? in 2015 (Fig. 1b). Even though the individual water 
footprints of the IWF, DWF and AWF increased, their percentages in the total amount of the total 
water footprint declined because of the very large increase in the CWF. 
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Fig.1 Composition percentages and amounts of the ecological footprint (a) and the water footprint (b) in Xinjiang 
from 2001 to 2015. IWF, industrial water footprint; DWF, domestic water footprint; AWF, animal-related water 
footprint; CWF, crop water footprint; TWF, total water footprint. The red line represents the variation trend. 


The ecological footprint varied significantly among prefectures and cities in Xinjiang (Fig. 2). 
Using 2015 as an example, Changji had the highest ecological footprint (31.18x10!° ghm’), 
followed by Aksu (16.89x10!° ghm?), Urumqi (16.70x10!° ghm?), Kashgar (15.83x10!° ghm?), 
Karamay (12.20x10!° ghm?), Hami (11.64x10!° ghm?) and Ili (11.32x10!° ghm?). Kizilsu had the 
smallest ecological footprint (1.10x10!° ghm?), as shown in Figure 2. 
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Fig. 2 Water and ecological footprints of the 14 prefectures and cities in Xinjiang in 2015. 1, Aksu Prefecture 
(abbreviated as Aksu); 2, Altay Prefecture (Altay); 3, Bayingol Mongolian Autonomous Prefecture (Bayingol); 4, 
Bortala Mongolian Autonomous Prefecture (Bortala); 5, Changji Hui Autonomous Prefecture (Changji); 6, Hami 
Prefecture (Hami); 7, Hotan Prefecture (Hotan); 8, Kashgar Prefecture (Kashgar); 9, Karamay City (Karamay); 10, 
Kizilsu Kirgiz Autonomous Prefecture (Kizilsu); 11, Tacheng Prefecture (Tacheng); 12, Turpan City (Turpan); 13, 
Urumqi City (Urumqi); 14, Ili Kazak Autonomous Prefecture (Ili). 


Figure 2 also illustrates the contrasting relationships between the water and ecological footprints 
of the 14 prefectures and cities. In some prefectures (Aksu, Altay, Bayingol, Bortala, Hotan, 
Kashgar, Kizilsu, Tacheng, Turpan and Ili), the ecological footprint was relatively small, but the 
water footprint was relatively large. In contrast, in Changji, Hami, Karamay and Urumqi, the 
ecological footprint was large, while the water footprint was relatively small. In Karamay and 
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Urumqi, the ecological footprints were relatively large but the water footprints were small, because 
the secondary and tertiary industries were relatively developed in these two cities, while other 
prefectures and cities were mainly dominated by agricultural production and thus large amounts of 
irrigation water were used, resulting in the large water footprints. 


4.2 Coupling relationship between the water and ecological footprints 


To examine the coupling relationship between the water and ecological footprints, we processed 
the two sequences as described in Section 3.3. As shown in Figure 3, both homogeneous indices 
for the water and ecological footprints showed increasing trends from 2001 to 2015, and the 
increase in the ecological footprint was larger than that of the water footprint. Specifically, the 
homogeneous indices of the water footprint increased from 0.73 in 2001 to 1.40 in 2015, with an 
increasing amplitude of 0.67 and obvious inter-annual variations. However, the homogeneous 
indices of the ecological footprint increased from 0.45 in 2001 to 1.71 in 2015, with an increasing 
amplitude 1.26, a level almost twice the increased amplitude of the water footprint. The coupling 
degree between the water and ecological footprints in Xinjiang during the study period remained 
extremely high, with an average coupling degree value of 0.99 and a standard error of 0.01. 
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Fig. 3 Homogeneous indices of the water and ecological footprints and coupling degree between the two 
footprints in Xinjiang from 2001 to 2015 


Spatially, the coupling degree between the water and ecological footprints in the 14 prefectures 
and cities varied significantly (Fig. 4). Except for Karamay and Urumqi, the other prefectures and 
cities exhibited a relatively high coupling degree (>0.80). Thus, in these prefectures and cities, the 
ecological footprint changed at the same pace as the water footprint, i.e., prefectures and cities that 
had a large (or small) ecological footprint also had a large (or small) water footprint. The coupling 
degree was relatively low (0.67) in Urumqi, thus the ecological footprint had a low coupling 
relationship with the water footprint. For Karamay, the coupling degree value of 0.43 indicated that 
the ecological footprint was not coupled with the water footprint. 

The low coupling degree values between the water and ecological footprints in Urumqi and 
Karamay indicated that the more developed a prefecture or city, the lower its water consumption. 
Agriculture in Urumqi and Karamay is underdevelopment; therefore, the CWF values were smaller 
than those in other prefectures and cities. In contrast, the urbanization and industrialization levels 
were relatively high in Urumqi and Karamay, resulting in large ecological footprints. 


4.3 Internal mechanisms between the water and ecological footprints 


The Granger causality test reflected the interactions and causal relationships between the water 
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Fig. 4 Coupling degree between the water and ecological footprints in the 14 prefectures and cities in Xinjiang 
in 2015. The numbers of the prefectures and cities are the same as in Figure 2. 


and ecological footprints (Table 2). Two different hypotheses (i.e., "the water footprint is not the 
Granger cause of the ecological footprint" and "the ecological footprint is not the Granger cause of 
the water footprint") and five lagging periods were tested in the Granger causality analysis shown 
in Table 1. According to the principles of hypothesis testing, when the significance reached the 0.05 
level (P<0.05), the null hypothesis was rejected. As shown in Table 1, the hypothesis "the water 
footprint is not the Granger cause of the ecological footprint" was rejected when the lagging periods 
were | and 2 years, thus leading to the conclusion that the increased water footprint was the Granger 
cause of the increased ecological footprint. This finding verified that an increase in the water 
footprint will cause an increase in the ecological footprint, but the effect emerges 1 and 2 years 
later. 


Table 2 Granger causality test of the water and ecological footprints 
Lagging : F Significance Coupling 
period Hypothesis statistic level decision 


The ecological footprint is not the Granger cause of the 


0 s 4.3 0.05 Accept 
water footprint 

0 The water footprint is not the Granger cause of the 08 0.36 Nipi 
ecological footprint 

1 The ecological footprint is not the Granger cause of the 34 0.05 Accept 
water footprint 

1 The water footprint is not the Granger cause of the 74 <0.01 Decline 
ecological footprint 

2 The ecological footprint is not the Granger cause of the 12 034 Avene 
water footprint 

2 The water footprint is not the Granger cause of the 52 001 Decline 
ecological footprint 

3 The ecological footprint is not the Granger cause of the 20 0.15 Accept 
water footprint 

3 The water footprint is not the Granger cause of the 20 0.16 Accept 
ecological footprint 

4 The ecological footprint is not the Granger cause of the 24 0.12 Accept 
water footprint 

4 The water footprint is not the Granger cause of the 12 0.37 Accept 


ecological footprint 
Note: 0, current year; 1, one year later; 2, two years later; 3, three years later; 4, four years later. 


The Granger causality analysis illustrated that in Xinjiang, the development of the social 
economy was dependent on the water consumption. However, this relationship did not always mean 
that the greater water consumption was associated with the higher level of social economy. For 
instance, Urumqi and Karamay were the two cities with the highest values of social economy in 
Xinjiang, but the water consumption and the water footprint were not very large, as compared with 
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other prefectures and cities. The modest water footprints occurred because these two cities had high 
levels of industrialization and mainly focused on the secondary and tertiary industries. In contrast, 
the social economies in other prefectures and cities in Xinjiang were still dependent on agriculture, 
which consumed more fresh water resources than industries. Thus, to decrease water consumption 
and preserve water resources in the arid environments in Xinjiang, it is extremely important to take 
measures to accomplish industrial transformation from the primary industry to the secondary and 
tertiary industries. 


5 Discussion 


As shown in Figure 2, prefectures and cities such as Changji, Hami, Karamay, Turpan and Urumqi 
occupied a relative high ecological footprint but a low water footprint (shown as the dot of the 
water footprint was in the column of the ecological footprint in Fig. 2), this was because these 
prefectures and cities were relatively more developed in the secondary and tertiary industries The 
secondary and tertiary industries have advantages of high economic benefit and low water 
consumption. For example, according to the Xinjiang Statistical Yearbook in 2016, the personal 
incomes (13.10x10* CNY) of Karamay was 13 times higher than that of Hotan (1.01x10* CNY), 
but the per capita water footprint in Hotan (1.10x10° m°) was more than three times of that in 
Karamay (0.3210? m3). Under this condition, the underdeveloped prefectures and cities intend to 
expand the scale of agriculture to bring more incomes. The agriculture consumes more fresh water 
resources, which will lead to environmental degradation; the poor environment of these prefectures 
and cities has low attraction to external investment, and less external investment enlarges the 
incomes gap between the low and high developed prefectures. This is a kind of Matthew effect, 
i.e., the rich getting richer and the poor getting poorer. Therefore, changing the structure of 
economy is the only way to accomplish harmonious and sustainable development of economy, 
environment and natural resources. 

The ecological footprint represents the land area required to sustain the current levels of resource 
consumption and the wastes generated by a population or individual. On the one hand, the 
ecological footprint indicates the influence of human behavior on the natural environment; on the 
other hand, it reflects the development of society and economy. How does the ecological footprint 
contribute to the local economy? Can development that is sustainable from an economic perspective 
be used to estimate the increase in the ecological footprint (ecological degradation)? To answer 
these questions, we determined the production value of the ecological footprint through dividing 
the regional GDP by the amount of the ecological footprint (Fig. 5). 
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Fig.5 Output value per unit ecological footprint in Xinjiang from 2001 to 2015 


Figure 5 shows that in Xinjiang, the output value per unit ecological footprint increases from 3.8 
to 6.4 CNY/ghm? from 2001 to 2015. This phenomenon illustrates that the output benefit of 
resource utilization increases with the development of social economy. However, the high resource 
usage efficiency cannot offset the effects of ecological damage and resource consumption resulting 
from the increase in the total ecological footprint. Previous research has shown that the expansion 
of the ecological footprint does not decrease the ecological carrying capacity, and the total carrying 
capacity is nearly stable (Xue et al., 2009). However, in the calculation of the ecological carrying 
capacity, only the area of land use was considered, resulting in the misunderstanding of the carrying 
capacity. For example, for the degenerated forests, if the area does not change, the calculated 
carrying capacity based on land use remains the same. However, it should be noted that the actual 
carrying capacity decreases. 

In this study, we further examined the analysis of homogeneous indices of the water and 
ecological footprint to better understand the relationship between the increases in the water and 
ecological footprints (Fig. 6). 
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Fig. 6 Relationships between the homogeneous indices of the water and ecological footprints (a) and between the 
water and biomass footprints (b) in Xinjiang from 2001 to 2015. The biomass footprint is the ecological footprint 
of renewable resources. 


Figure 6a illustrates that the increase in the ecological footprint is larger than that in the water 
footprint. In Figure 6b, the biomass footprint refers to the consumption of renewable resources 
(calculated as the ecological footprint minus the footprint of fossil fuels as stated in Section 3.1), 
and it reflects the pressure of bio-productive requirement on the environment (Feng, 2011). The 
linear regressions between the homogeneous indices of the water and ecological footprints and 
between the water and biomass footprints show that as the water footprint increases, the ecological 
and biomass footprints rise up correspondingly, and the slopes of the regressions are 2.64 and 1.91, 
respectively. In other words, the regression equations indicate that when the water footprint 
increased by 1 unit, the ecological and biomass footprints increased by almost 3 and 2 units, 
respectively. That is to say, l-unit increase in the water footprint would result in approximately 2- 
unit increase in the biomass footprint or 3-unit increase in the ecological footprint. Thus, we 
proposed a hypothesis that an increase in the water resource usage in the artificial oasis can result 
in a twice to triple of the ecological degradation in the natural oasis. The logical conclusion from 
this finding is that the formation and expansion of 1 unit of the artificial oasis will cause the 
degradation of 2—3 units of the natural oasis. 

In order to validate this hypothesis, we took the Tarim River Basin as an example and try to 
contrast the water consumption in both natural and artificial oases. The water consumption and the 
area of the natural oasis were obtained from Ling et al. (2014). We calculated the water 
consumption in the artificial oasis by statistic data from the Xinjiang Water Resources Bulletin 
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(Water Resources Bureau of Xinjiang Uygur Autonomous Region, 2001-2015), and derived the 
areas of the artificial oasis in Tarim River Basin from the aforementioned land use in 2010 
according to Fan et al. (2002) and Wang et al. (2018). 

Ling et al. (2014) has calculated the ecological water consumption of the natural oasis in the 
mainstream of Tarim River in the year of 2010 based on the remote sensing data. According to this 
previous result, we extracted the area of natural vegetation and the amount of water consumption 
of the natural oasis in different sections and regions in the mainstream of Tarim River, and then 
calculated the water consumption rate, as shown in Table 3. 


Table 3 Area and water consumption of the natural oasis in the mainstream of Tarim River 


Section 1: Alar—Xinquman Section 2: Xinquman—Yingbazha 
Region anes Water. k e prea ot Nater Pii a 
vegetation consumption vegetation consumption 
(x10? km?) (x108 m’) at Ps aes (x10? km?) (<108 m3) a Fe 
Nae 2.27 3.84 1.69 3.11 5.75 1.85 
aie 1.66 1.86 1.12 1.53 1.83 1.20 
Total 3.93 5.70 1.45 4.64 7.58 1.63 
Section 3: Yingbazha—Wusiman Section 4: Wusiman—Qiala 
Region Aten o1 Nater a Ateno Laat ene 
vegetation consumption vegetation consumption 
(x10? km?) (x108 m°) a P n (x10 km?) (x108 m°) i 6 m/km?) 
neg 2.86 3.83 1.34 2.12 3.80 1.79 
iat 0.94 1.10 1.16 0.30 0.30 1.01 
Total 3.80 4.93 1.30 2.42 4.10 1.70 
Section 5: Qiala—Taitema Lake Tarim River 
Region Atea OF Water M ATTOL Water ee 
vegetation consumption vegetation consumption 
(«10* km?) CIm) caom) | C10 km?) («108 m°) (10? men?) 
pare 0.93 1.33 1.43 11.29 18.56 1.64 
sie 0.57 1.21 2.14 5.00 6.30 1.26 
Total 1.50 2.54 1.70 16.29 24.86 1.53 


According to Table 3, the area of natural vegetation in the mainstream of Tarim River was 
16.29x10° km?, and the total amount of water consumption was 24.86x108 m?. Due to the 
differences in the vigor and composition of natural vegetation among the sections (Ling et al., 
2014), the water consumption rate varied from 1.01*10° to 2.14 10° m?/km?. 

In order to contrast with the water consumption in the artificial oasis, we analyzed the area and 
the amount of water consumption of the artificial oasis in the five prefectures in Tarim River Basin. 
The results are also shown in Table 4. Kizilsu had the smallest artificial oasis area (0.8010? km?) 
but the highest water consumption rate (6.85105 m?/km7). This was due to that Kizilsu was popular 
in livestock husbandry and the water consumption of this region accounted for 54% of the national 
water consumption. But in other prefectures, the irrigation water was the main source of water 
usage, varying from 55% to 70% of the regional total water consumption, and the water 
consumption rate varied from 2.32105 to 4.40x10° m?/km?. To sum up, the ranges of the water 
consumption rate in the artificial and natural oases were 2.32x10°-6.85x10° and 1.01x10°— 
2.14105 m?/km?, respectively (Tables 3 and 4). 

As mentioned above, the water consumption rate in the artificial oasis was 2—3 folds of that in 
the natural oasis. This result supports the hypothesis that the formation and expansion of | unit of 


the artificial oasis will cause the degradation of 2—3 units of the natural oasis. It should be 
emphasized that the existence of the natural oasis has particular importance. First, the natural oasis 
is home to numerous species of flora and fauna, and many of which are rare; thus, to protect the 
biodiversity, it is of importance to protect the natural oasis. Second, the natural oasis is a defensive 
‘shield' for humans who live there. The natural oasis surrounds the artificial oasis, thereby 
minimizing the exposure of the artificial oasis to the desert and providing a comfortable living 
environment. Thus, the protection of the natural oasis is to protect the living conditions of humans. 


Table 4 Area and water consumption of the artificial oasis in the five prefectures in Tarim River Basin 


Prefecture Area : Water consumption Water consumption rate 
(x10? km?) («108 m°) (x105 m3/km’) 

Bayingol 12.33 28.66 2.32 
Hotan 12.74 33.95 2.67 
Aksu 16.19 71.23 4.40 
Kashgar 18.44 13399 3.99 
Kizilsu 0.80 5.46 6.85 
Total 60.50 212.85 3.52 


According to the analysis of this research, the increase in the water consumption was the main 
cause of ecological deconstruction in Xinjiang. However, the cases of Urumqi and Karamay 
illustrate that the development of a city does not inevitably cause the increase of water 
consumption; that is, economic development does not always have a direct relationship with 
ecological degradation. Therefore, for the artificial oasis, the structure of economy should be 
transformed by enhancing technology to solve the problem of short-term overuse of resources. 

The ecological footprint assessment is an efficient tool to evaluate the sustainability of social- 
economy and environment. While there still be some controversies and criticism on this method. 
In the progress of the ecological footprint calculation, equilibrium factors and yield factors are 
uniformed and offered by Wackernagel (1998) as global average value. However, the characteristics 
of land resources vary significantly among different regions, which indicates that the use of uniform 
values cannot really reflect the actual situation of the ecological productive land in a region, and it 
will mislead the sustainable use of resources (Li et al., 2011; Dai et al., 2015; Somogyi, 2016). In 
addition, the calculation of the ecological footprint was based on an ideal hypothesis that land 
resources merely have a single function, while in reality, land resources sometimes have more than 
one utility. For instance, in some regions, crops could be harvested for twice per year or three times 
every two years, but this situation is not considered in the ecological footprint model, resulting in 
errors in the calculation results. 

The assessment of the water footprint has advantage of evaluating the spatial distribution. 
However, most of the previous researches cannot accurately reflect the changes of the global water 
footprint due to the lack of detailed spatial geographic features (Mekonnen and Hoekstra, 2014). 
Further, it is hardly to apply the water resource management accurately and effectively (Perry, 
2014). Therefore, we suggest that researchers should considered the time and space information in 
the water footprint assessment in the future, in order to accurately evaluate and compare the 
consumptions of water resources by human activities. The global virtual water trade and water 
footprint in the production process should be accounted for separately, because of the different 
opportunity costs and impacts of the green, blue and grey water footprints (Hoekstra and 
Chapagain, 2008). In addition, the CWF merely accounts the water consumption and pollution, and 
it cannot reflect the impact of water consumption on the environment directly (Gao et al., 2017). 
Therefore, it has a great significance to reveal the deeply interaction between the water and 
ecological footprints. 


6 Conclusions 


In this study, the variations in the water and ecological footprints in the 14 prefectures and cities in 
Xinjiang of China from 2001 to 2015 were evaluated. The integrated relationship between the water 
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and ecological footprints was analyzed with the Granger causality test. The relationship between 
the increases in the water and ecological footprints was quantified. 

Both the water and ecological footprints increased significantly from 2001 to 2015. Specifically, 
the water footprint increased from 293.12x108 to 558.47x108 m? and the ecological footprint 
increased from 38.80x10!° to 146.94x10!° ghm?. In general, the ecological footprint was highly 
coupled with the water footprint both at temporal and spatial scales, because most of the prefectures 
and cities in Xinjiang were dominated by agriculture. However, in the highly developed cities 
(Urumqi and Karamay,) the water and ecological footprints were relatively uncoupled. The increase 
in the ecological footprint was stimulated by the water footprint, but the effect lagged by 1 and 2 
years. The increase in the ecological footprint was 2-3 folds that of the water footprint, and the 
water consumption rate in the artificial oasis was 2—3 folds that in the natural oasis; thus, we deduce 
that the formation and expansion of | unit of the artificial oasis would lead to the degradation of 
2-3 units of the natural oasis. 
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